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■ Abstract 

The pion and nucleon mass differences generate a very pronounced cusp in the pho- 
toproduction reaction of a single pion on the nucleon. A nonrelativistic effective field 
theory to describe this reaction is constructed. The approach is rigorous in the sense 
that it is an effective field theory with a consistent power counting scheme. Expres- 
sions for the S- and P-wave multipole amplitudes for all four reaction channels at 
two loops are given. 
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CO 
ON 

1. Isospin breaking effects have recently received a considerable amount of 
attention. This was triggered by the observation of a cusp-like structure in 
the decay rate of K — > 3tt decays [1], which originates from isospin break- 
ing. Since the strength of this cusp is intimately related to the 7T7t scattering 
5_i ■ lengths, another possibility to measure the 7T7r-scattering lengths - beside K e 4 



decays [2, 3] and the lifetime of the pionium atom [4] - was established. It 
is amusing to note that the same effect that is responsible for the cusp in 
K — > 3tt decays also leads to a correction to the 7T7r-scattering phase shifts 
from K ei decays [5]. 

The photoproduction reaction of neutral pions on the proton is a reaction 
where isospin breaking corrections have not yet been calculated in a fully sys- 
tematic approach. On the other hand, it has been known for a long time that 
this reaction shows a very strong effect due to isospin breaking: the electric 
multipole E 0+ exhibits an exceptionally strong cusp at the ir + n threshold (see 
for instance Ref. [6]). In much the same way as in K — > 3ir decays, a nonrela- 
tivistic theory - adapted in this letter to the photoproduction reaction of pions 
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on nucleons - provides a rigorous framework which yields quantum field the- 
oretical matrix elements with the effective range parameters of pion-nucleon 
scattering and the threshold parameters of the photoproduction reaction as 
free parameters. Furthermore, the framework allows one to include radiative 
corrections in a standard manner at a later stage. 

The cusp in neutral pion photoproduction has been studied before. Ref. [7] 
introduces a two-parameter model, which captures the most important leading 
effect of the cusp. In Ref. [8], a coupled channel S- matrix approach is used to 
investigate the cusp structure. 

In a first step, the nonrelativistic Lagrangian is written down and the power 
counting is discussed. Then, the coupling constants are matched to the perti- 
nent threshold parameters. The calculation of the multipole amplitudes up to 
two-loop order is then straightforward. The representation of the amplitudes 
correctly reproduces the analytical structure in the low energy region. As an 
application, the relation of the phase of the electric multipole E 0+ to the phase 
of the S-wave of rr°p — > n°p scattering is discussed. 

2. Some basic relations and definitions used in the analysis of pion photo- 
production are collected. We calculate the matrix element for the process 
N(pi) + j(k) — » N(p 2 ) + 7r 0,± (g) for all four channels at leading order in the 
electromagnetic coupling e, 

(p 2 , q out|pi, k in) = - ? (2vr) 4 (5 (4) (P/ - P) u(p 2 , t')e^u( Pl ,t) , (1) 

where N denotes either a proton or a neutron, P; and Pf denote the total 
four momentum in the initial and in the final state, respectively and & stands 
for the polarization vector of the photon. In the following, the four reaction 
channels will be abbreviated as 

7P — > pn° : (pO) , 7J9 — > nn + : (n+) , — > nn° : (nO) , — > pn~ : (p— ) . 

To analyze the photoproduction reaction of pions, electric and magnetic mul- 
tipoles are usually introduced. To this end, the amplitude is written in terms 
of two component spinors £ t and Pauli matrices r k [9] , 

A^ = 87r v ^^9 r 6, 
J = ir • e ?i + t ■ q r • (k x e) ? 2 + «t ■ k q ■ e J 3 + ir ■ qq ■ e ? 4 . (2) 

The hat denotes unit vectors. The 3^ are decomposed into electric and mag- 
netic multipoles with the help of derivatives of the Legendre polynomials 

Pi{z)[% 
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?i = £ [ZM,+ + £ l+ ]P/ +1 Cz) + [(/ + 1)M,_ + E^)PU{z) , 



1=0 



? 2 = Y,[(l + l)Mi+ + lMi-)P!{z) , 



1=1 
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1=1 



?4 = E [ M *+ - ^+ - M_ - (*) • (3) 

1=1 

The discussion is restrained to the center of mass frame in the rest of the 
article. 

3. To describe the behavior of the multipoles close to threshold - where the 
energy of the produced pion and of the nucleon are small - a nonrelativistic 
calculation is justified. Furthermore, it offers the advantage that all the masses 
can be set to their physical value. Therefore, all the poles and branch points 
appear at the correct place in the Mandelstam plane. Moreover, the interaction 
of the nucleon and the pion is described by effective range parameters, which 
allows one to directly access the pion-nucleon scattering lengths. 

The covariant formulation of nonrelativistic field theories introduced in 
Refs. [10-12] is used here since it incorporates the correct relativistic disper- 
sion law for the particles. The nonrelativistic proton, neutron and pion fields 
are denoted by ■0, x an d ^k, respectively. The kinetic part of the Lagrangian 
after minimal substitution takes the form (see Ref. [13]) 

C kin = W«4AW ± 7T ± - i(AW ± 7r±)t7r ± - 27riW 2 ± 7r ± 
± ^ 

+ 2 X ] W n {id t - W n ) X + 27rlW (id t - W )tt , (4) 

with 



Wo = JM* - A , W n = J ml - A , D t n ± = (d t t teA )ir± 



D t ip = (d t - leAo)^ , W± = yJM* — D 2 , W p = A /m 2 - D 



Y2 



Dvr± = (V ± ieA)TT± , = (V + ieA)ijj . (5) 

Note that since the photon is treated as an external field, its kinetic term is 
absent. 

4. Consider the reaction^ p(pi) + 'j(k) — > p(p2) + vr°(g). Close to threshold, 
the momenta of the incoming proton and photon are of the order of the pion 



The arguments remain of course the same for the other three channel. 
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mass whereas the outgoing particles have very small momenta. Therefore, 
we count momenta of the outgoing pion and the outgoing proton as a small 
quantity of 0(e) and the momenta of the incoming proton and of the photon 
as 0(1). All the masses are counted as 0(1). The mass differences of the 
charged and neutral pion, A^ = — M^ and of the proton and the neutron, 
A at = m^—rrip are counted as 0(e 2 ). We refrain from counting the pion mass as 
a small quantity with respect to the nucleon mass, since considering a particle 
as light and nonetheless describing it nonrelativistically seems inconsistent 
(see also Ref. [14]). 

At first sight, this counting scheme seems to lead to infinitely many terms 
already in the leading order p + 7 — * p + n° Lagrangian £ 7 because derivatives 
on the incoming fields are not suppressed in e. However, since the the modulus 
of the momentum of the incoming particles, |k|, can be expanded in the small 
momentum |q|, 

| k | = W ", fc = ^l±», * = £ + y=*Z, (6) 

n 2 1 + V 4M n o(l + y) m p 

one obtains a valid power counting scheme: Consider the Feynman rule in 
momentum space of an operator of order 0(e°) with a given arbitrary number 
of derivatives. This expression can be expanded in a sum of one term of order 
0(e°) without any momenta of the incoming particles present and subsequent 
higher order terms. Doing this for every operator of 0(e°), all the resulting 
leading order terms without any momentum dependence can be described by 
one operator of order 0(e°) in the interaction Lagrangian. The same proce- 
dure leads to finite numbers of operators at any given higher order in e. The 
derivatives on the incoming fields are only needed to generate unit vectors 
in the direction of the incoming photon. This shows that the nonrelativistic 
theory is not capable of predicting the dependence on |k| even at threshold. 
An additional generic parameter a is introduced to count the pion-nucleon 
scattering vertices. Every pion-nucleon interaction vertex counts as a quantity 
of order 0(a) since the coupling constants are proportional to the pion-nucleon 
scattering threshold parameters, which are small. The perturbative expansion 
is therefore a combined expansion in e and a. 

5. The Lagrangian needed for the calculation of the amplitudes for pion pho- 
toproduction reads C = Ckin + £7 + £wN , where Chin denotes the kinetic part, 
£ 7 incorporates the interaction with the photon field and and C n N describes 
the pion-nucleon sector. 

In the pion nucleon sector, the leading terms of the Lagrangian have been 
given before in Ref. [15]. First, some notation is introduced in order to write 
the Lagrangian in a compact form. For every channel n, we collect the charges 
of the outgoing and the incoming pions in the variables v and w, (n;v,w): 
(0; 0, 0), (1; 0, +), (2; +, +), (3; 0, 0), (4; -, 0), (5; -, -), thereby assigning 
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unique values to the variables v and w once n is given. The Lagrangian reads 
(, 



[Tl A} T {2 , 3} ) \ x )' 



(7) 



neC 

with the abbreviation / A g = fAg + (Af)g. 

For £ 7 , the photon is treated as an external vector field A which is odd under 
parity and time-reversal transformations. Gauge invariance requires that it can 
only appear in covariant derivatives and through the Maxwell equations in the 
electric and magnetic field E = — VA° — A and B = V x A. The interaction 
Lagrangian can be read off from the nonrelativistically reduced and multipole 
expanded matrix element, Eq. (2) and the known threshold behavior of the 
multipoles, E t ± : M t ± ~ \q\ l . As already mentioned, the effective theory fails 
to predict the dependence on |k| even at threshold. Factors of |k| can always 
be obtained by a redefinition of the coupling constants of the Lagrangian. 
Of course, all the terms have to be invariant under space rotations, parity 
and time reversal transformations. This leads to the following expression for 
the Lagrangian C 1 . The upper index on the coupling constants counts the 
number of derivatives on the external vector field and is introduced for later 
convenience. 

= -iG® ^r k i) V ] E k V%d + iGP ^r m r l i) B l X7 m 4 



-iGf ^r j ^V 3 E k v k 4 
4 2 ) = -iGf^T k ^V jl E k V jl 7xl - iG < i ) ^r k ijE k A^ 



7 

"G^ -yj- I ' -iy V " Hi V /Iq 

t •^yfl) ,i h 1 7-ifc A 

^0 

+ iG i e ) ^T m T l ^V n B l V mn 4 - iG^ ^ ] T^V jl E k V kl 4 

- iG < i ) ^T j ^E k V jk 4 

£(3) = -iG < i ) ^T k i)V ] E k AV ] iil - iG ( §^T k i)V lmn E k V lmn nl 
+ iG ( ^^T m T l ^B l AV rn 4 + iG ( $^T m T l ^V in B l V min 4 

- iG { S^r^V j E k AV k 4 - iG { $^T j 4>V jlm E k V klm 4 

- iG { S^r^V l E k V jkl 4 (8) 

Here, the notation V n * 2 '"* fc = V* 1 V* 2 • • • V ifc is used. Since the structure of the 
Lagrangian for the remaining channels stays the same, one only has to replace 
the coupling constants and the field operators, 

(n+) : {^,4,G\ n) } - {x*A,HV>}, (nO) : ty, ^, Gf} - {x,x\L^}, 
(p-):{1>AM n) }^{x,*UKl n) }. (9) 
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Fig. 1. The ratio R plotted as a function of the photon energy in the lab frame. The 
dashed line shows the nir + threshold. 

The full interaction Lagrangian £ 7 is then given by adding the £W of all four 
channels. 

Here, a remark about the structure of the interaction terms seems in order. 
Since the interaction Lagrangian contains by construction no nNN couplings, 
there is no nucleon pole diagram. Therefore one might be worried that an 
expansion in powers of the momentum q does not converge sufficiently fast 
even in a vicinity of the threshold. However, we check explicitly that the 
nucleon pole in the s channel can be approximated by a polynomial in q. 
The tree level expressions of E 0+ in the channel (pO) (taken from Ref. [7]) 
expanded up to and including terms of order q 2 yields a good approximation 
of the full tree level result. In Fig. 1, the ratio R = (i?o+ P ~~ E 0+ )/E 0+ is 
plotted as a function of the photon energy in the rest frame of the proton. In 
view of this result, we also conclude that there is no reason to be concerned 
about singularities from resonances like the p in the t channel or the A(1232) 
in the s channel, which are almost as close to threshold. 

6. In the pion-nucleon sector, the coupling constants of the nonrelativistic 

(k) 

Lagrangian, Cj and D\ can be expressed in terms of pion-nucleon scattering 
lengths of the S'-wave and P-wave, a 0+ and a±± and effective range parameters 
b 0+ , respectively. Adopting the notation of Ref. [16], in the isospin limit, the 
isospin decomposition of the nN scattering amplitudes reads 

2p7r°— ►jwr -^rar — m7r° T , T^ )7r o_ >n7r + T n7T Q^p w - \PlT , 

-^~?i7r+^n7r+ -^~p7r~ — *pn~ "F T . (10) 

Defining Af = 47i(m p + M n ), one finds 

C = 2 A/"4+ , Ci = 2v / 2A/a 0+ , C 2 = 2 A/"(4 + + % + ) , 

C3 = Co , C4 — C1, C5 = C2 ■ (11) 
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The matching conditions for the are given in a generic form only. The 
isospin index of the threshold parameters can be inferred from Eq. (10)[H 

= 2M(2a 1+ + fll _) , A (2) = -M { + bo+) , 



2m p Mv 



Df ] = 2Af (ai_ - a 1+ ) . (12) 

Here, higher order terms in the threshold parameters have been dropped. The 
corrections to these relations which appear due to isospin breaking have to be 
calculated within the underlying relativistic theory. For the C», they can be 
found in Refs. [17-19]. Note that the second line in Eq. (11) is only true in 
the isospin limit. 

The constants G[ n \ H\ n \ i*Q and on the other hand are related to the 
threshold parameters of the electric and magnetic multipoles of the pertinent 
channel. In the isospin limit, the expansion of the real part of the multipole 
Xi± close to threshold is written in the form 

oo 

ReX l± (s) = J2X l± , 2k \q.\ l+2k , (13) 

fc=0 

which defines the threshold parameters Xi±^k- In the following, the relations 
of the coupling constants G^ to these threshold parameters is given at lead- 
ing order in the pion nucleon threshold parameters. Since the nonrelativistic 
theory is not suited for the study of the dependence of the multipoles on |k|, 
in this analysis, all vectors k are turned into unit vectors by the pertinent 
redefinition of the coupling constants, 

Gf } = Af k n Gi , Af = Aix{m p + M w o) . (14) 

Note that the higher order corrections due to Eq. (6) are taken care of in the 
matching relations. Again, these relations pick up isospin breaking corrections 
which have to be evaluated in the underlying relativistic theory. 

Only the matching equations for the couplings of the Lagrangians and 
Cy are indicated in the main text, the remaining relations are relegated to 
appendix A. To ease notation, X i± = Aj ±j o is used. 

G = 2E 0+ , G 1 = 6(E +1 + M +1 ), 

G 2 = -2(M_! + 2M +1 ) , G 3 = 6(E 1+ - M 1+ ) . (15) 

For the coupling constants Hi, Ki and Lj, the algebraic form of the relations 
is identical. However, the multipoles of the pertinent channels appear and the 
masses in Eq. (14) have to be adjusted. 



Note that we use the Condon-Shortley phase convention. 
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Upon the inclusion of dynamical photons, the nonrelativistic coupling con- 
stants pick up an imaginary part due to 'yN intermediate states (see Ref. [14] 
for a discussion of this issue). For the coupling constants Cj, the imaginary 
part can be obtained from Ref. [19]. The imaginary part of the coupling con- 
stants Go and H in the photoproduction Lagrangian £ 7 follows from the 
imaginary parts of E 0+ at threshold. A calculation at leading order in chiral 
perturbation theory [20] yields 



Im£ 0+ - e ' 9V{v + 2) 



647r 2 F(l + ?/)- . 
for the channel (pO) and 
e 3 g(y + 2) 



(2/ + l) 3 ln(l + 2/) 



(16) 



lmE 0+ = 



V2 327r 2 F(l + y)^ 



-(2y + 5)(y 2 + 2)-(y + iy\n(l + y) 



(17) 



for the channel (n+). Here, F and g denote the pion decay constant and the ax- 
ial coupling, both in the chiral limit. A comparison with the real part at thresh- 
old - which is estimated taking the experimental results from Refs. [21,22] 
shows that the imaginary part is indeed of the generic order of electromagnetic 
corrections, 

Im Gn „/ o N Im Ho „ . . 

r^- < 10 >• r^t°< 10 >■ (18 » 

In the following, we assume that the coupling constants are real. 

7. In the following, we provide the expressions for the electric and magnetic 
multipoles Ei + for I = 0, 1 and Mj± for I — 1 for the channel (pO). The result 
is written in the form 

X l!± (s) = Xfr(s) + XS?*(8) + Xf^(s) ■■■ (19) 

where s = (pi + k) 2 and the ellipsis denote higher order terms in the expansion 
in e and a. The results for the other channels can be recovered by a simple 
replacement of the coupling constants which will be given later. Writing 

xt^(s) = x} ± a l + xf ±t2 d 2+l + ■■■. (20) 

one finds 

Eq + = Go , 3Eq + 2 = G4 — 3G5 + Gq — Gs , 

6M' + = G\ — G3 , M[ + 2 = ~^Gq + JqGio + J5G12 + |Ci3 — , 

3M*_ = G3 — G\ — 3G2 , Ml_ 2 = |Gg — |Gio + Gu — |Ci2 

— |Cl3 + J5G14 , 

&E\ + = Gi + G3 , ^.+,2 = ~~ \G§ + jqGio + j$Gi2 — \Giz 

+ ^Gi4 - ^Gi 5 . (21) 
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Fig. 2. One- and two loop topologies needed to calculate the amplitude. The dou- 
ble line generically denotes a nucleon, the dashed line a pion and the wiggly line 
indicates the external electromagnetic field. 

One observes that D-waves appear naturally at order e 2 in this framework 
(see also Ref. [23]). 

8. All the one-loop contributions are proportional to the basic integral 



i{2n) D 2u a {\)2u b {V - 1) (u a (l) - Z )MP -l)-Po + k) 



w±(p) = yJM* + P 2 , 
In the limit D — > 4, 



p2 



m 2 + p 2 



p2 



{22) 



J ab {P 2 ) 



16ns 



sj(s - (m a + M wb ) 2 ){s - {m a - M n t) 2 ) , (23) 



which is a quantity of order e. The one-loop result for channel (c) up to and 
including order 0(ae 4 ) reads 



\ 



|q| 



l_^-ilLoop ( 



1 + 



( p(c) p(c) \ 

-'ll ^12 

p( c ) p(c) 

^21 -*22 

p( c ) p( c ) 

31 -<32 

. pO) p(c) . 

\ -<41 -"42 / 




(24) 



where m a , M^i, denote the masses of the final state of the pertinent channel 
and m c , M n d stands for the masses of the intermediate state that differ from 
the final state masses. The elements are functions of the pion momentum 
q and the coupling constants of the Lagrangian. For the channel (pO), one 
finds 



p(po) 
*ii 

P (po) 

18P 2 ( ? 0) 
18P 2 ( f } 



i(2)/ 



,(2) 



,(2) 



(3) 



h (s, m c , M^) 



\G\ — G 3 ) , 
D? - D?) {Hi 
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9P 3 ( ? 0) = q 2 (d™ + 2D^) (G s - Gx - 3G 2 ) , 
gp(f ) = tf( s> mc , M wd ) (D? + 2D? >) (ff 3 - H x - 3H 2 ) , 
18P 4 r = q 2 - D®) ( Gl + G 3 ) , 

18pW = h 2 (s, m c , M nd ) (D? - D?) (H, + H 3 ) , (25) 

where Eq+* 2 denotes the pertinent coefficient of the tree level result of channel 
(c), see Eq. (20), and h 2 (s,m c , M nd ) is given by 

/i (s, m c , M nd ) = — , (26) 

which is a quantity of order e 2 . Eq. (24) and (25) clearly show the advantage 
of the nonrelativistic description: The strength of the cusp in the channel (pO) 
is parameterized in terms of the coupling constant C\ and the ratio H /G . 

9. The two-loop corrections all have the topology shown in Fig. 2 and can 
therefore be cast into the form 

(%%) (;;:). ™ 

where the for the channel (pO) read 
T[f > = G 2 G + Cg^SgV - 4G Go4V , 

Ti ( 2 0) = ±C? Go + \C, d ff + ^ < 2 ' V - Ci ff 2?SV - d G 

- l -C Q H D {2) q 2 + l -C Q G x h 2 (s, m c , M n d) 

- d Go £>J 2) /i 2 (s, m c , M nd ) - ^Go # D {2) h 2 (s, m c , M nd ) , 

T 2 { f = C&Ho - C 2 H,D {2 \ 2 + G 1 G 2 ^ 2 ) ' i /i 2 (s, m c , M wd ) 

- C 2 H D {2) h 2 (s, m c , M nd ) - 2C 1 H Q D { 2 ) h 2 (s, m c , M nd ) . (28) 

Up to and including order 0(a 2 e 5 ), the two-loop corrections are independent of 
the scattering angle cos and only contribute to the amplitude jFi. Therefore, 
at the order considered here, the electric multipole Eq + is the only quantity 
which receives two-loop corrections. 



10. The result for the other channels are obtained from the tree level result 

t 0) and T<f 



of channel (pO) and the coefficients Pjf® and Tjf® by the replacements 



(n+) 
(nO) 

(P-) 



{Gi, Hi, G , G 2 } — > {i^j, Gj, G 2 , G } , 

{Ci, ifj, Go, Gi, G 2 } — > {Li, Ki, G3, G4, G5} , 

{Gi, Hi, G , , Gi, G 2 } — > {.fQ, Lj, G 5 , G 4 , G 3 } . (29) 
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The replacement indicated for the C x has to be done also for the corresponding 

11. In an isospin symmetric world, the phase of the multipole E 0+ is directly 
related to the phase shift of the >S-wave of pion-nucleon scattering by virtue of 
the Fermi- Watson theorem [24]. In Ref. [8], it is shown with a coupled channel 
S'-matrix approach that in the channel (pO), to leading order in e, below the 
n + n threshold, the phase of the S wave of rr°p — > n°p scattering is equal to 
the phase of E 0+ , 

tano p7r o^ p7r o = tan— — - — = tano£ 0+ . (30) 
Re Eq + 

The framework developed here allows one to test this statement order by order 
in the perturbative expansion. To this end, the phase of E 0+ below the second 
threshold is calculated up to and including 0(a 2 e 4 ), 

tan5 Bo+ = C ImJpo + C\ J n+ lmJ p0 - 2L> 2) ImJ p0 q 2 (31) 
- 2Ci£>} 2) J„+ImJ p0 q 2 - 2dDf ] J n+ lmJ p0 h 2 (s, m n , M w ) + • • • . 

Calculating n°p — > 7i°p scattering to the same order with the Lagrangian given 
in Eq. (7), one finds that the phase of the S- wave below the second threshold 
is indeed equal to Eq. (31). However, the main object of interest is the phase 
of the S'-wave of 7i°p — > n°p scattering in the isospin symmetry limit, 

C 2 a' 2 

tan 8 pn o^ pn o = C lmJ p0 + -f ImJ n+ ~ a^ + q + H , (32) 

o a 0+ 

which does not agree with 8 p7T o^ p7T o in the presence of isospin violations already 
at leading order. 

12. In this letter, the photoproduction reaction of pions on the nucleon is 
studied using a nonrelativistic framework. The electric and magnetic multi- 
poles Ei + for I — 0, 1 and Mi± are calculated in a systematic double expansion 
in the final state pion- and nucleon momenta (counted as a small quantity of 
order e) and the threshold parameters of nN scattering (denoted by a). Ex- 
plicit representations for the multipoles up to and including e 3 , e 4 a, e 4 a 2 are 
provided. The representation is valid in the low energy region, at least up to 
a photon energy in the lab frame of E 1 = 165 MeV. It accurately describes 
the cusp structure and allows one to determine the pion-nucleon threshold 
parameters from experimental data. 

The relation of the phase of the electric multipole E Q+ in the (pO) channel to 
the phase of the >S-wave of 7r°p — > 7r°p scattering is discussed in the presence 
of isospin violation. A relation found in earlier work [8] is confirmed. We 
stress that the relation does not allow one to obtain the phase of n°p — > n°p 
scattering in the isospin limit. 
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A Matching relations 

The matching relations of the nonrelativistic couplings to the threshold pa- 
rameters defined in Eq. (13) read 

G 4 = 15{E 2+ + 2M 2+ ) , 

G 5 = 2(E 0+ - 2£ 2 _ - kE 0+ - 2E 0+)2 + 3{E 2+ - 2M 2 _ + 2M 2+ ) , 

G e = -6(3M 2+ + 2M 2 „) , 

G 7 = 30(E 2+ - M 2+ ) , 

G 8 = 6(£ 2 _ - M 2+ + M 2 _ + E 2+ ) , 

G 9 = 2d( - 6E 3 _ - 3kE 1+ - 6E 1+ , 2 + IQE 3+ - 24M 3 _ - 3kM 1+ 

- 6M 1+i2 + 45M 3 + , 
G w = 35(E 3+ + 3M 3+ ) , 

Gii = G 2 C, + kMi_ + 2M X _ 2 - 9M 3 _ + 2kM 1+ + AM 1+a - 12M 3+ , 
Gi 2 = -15(3M 3 _+4M 3+ ), 

G 13 = 2G 3 C - QE 3 _ - 3kE 1+ - QE 1+)2 + 15E 3+ - 6M 3 _ + 3kM 1+ 

+ 6Mi +)2 - 15M 3+ , 
G u = 105(E 3+ - M 3+ ) , 

G15 = 30(M 3 _ + £ 3 _ + E 3+ - M 3+ ) , (A.l) 

With K = M^~ P and C = V - 



References 

[1] J. R. Batley et al. [N A48/2 Collaboration], Phys. Lett. B 633 (2006) 173 
[arXiv:hep-ex/0511056] . 

[2] S. Pislak et al. [BNL-E865 Collaboration], Phys. Rev. Lett. 87 (2001) 221801 
|arXiv:hep-e x/0106071 1 . 

[3] J. R. Batley et al. [NA48/2 Collaboration], Eur. Phys. J. C 54 (2008) 411. 

[4] B. Adeva et al. [DIRAC Collaboration], Phys. Lett. B 619 (2005) 50 
|arXiv:hep-ex/0504044] . 



12 



[5] G. Colangelo, J. Gasser and A. Rusetsky, Eur. Phys. J. C 59 (2009) 777 
[arXiv:0811.0775l [hep-ph]]. 

[6] A. M. Bernstein, M. W. Ahmed, S. Stave, Y. K. Wu and H. R. Weller, 
larXiv:0902.3650l [nucl-ex]. 

[7] V. Bernard, N. Ka iser and U. G. Meissner, Z. Phys. C 70 (1996) 483 
[arXiv:hep-ph/9411287j. 



[8] A. M. Bernstein, Phys. Lett. B 442 (1998) 20 |arXiv:hep-ph/ 9810376 1 . 

[9] G. F. Chew, M. L. Goldberger, F. E. Low and Y. Nambu, Phys. Rev. 106 
(1957) 1345. 



[10! 

[11 
[12 

[13! 
[14 
[15! 

[16! 
[17! 

[18 

[19 

[20! 
[21 
[22 
[23 



G. Colangelo, J. Gas ser, B. Kubis and A. Rusetsky, Phys. Lett. B 638 (2006) 
187 |arXiv:hep-ph/0604084| . 

M. Bissegger, A. Fuhrer, J. Gasser, B. Kubis and A. Rusetsky, Phys. Lett. B 
659 (2008) 576 [ arXiv:0710.4456l [hep-ph]]. 

B. Kubis and S. P. Schneider, Eur. Phys. J. C 62 (2009) 511 |arXiv:0904.1320l 
[hep-ph]]. 

M. Bissegger, A. Fuhrer, J. Gasser, B. Kubis and A. Rusetsky, Nucl. Phys. B 
806 (2009) 178 [ arXiv:0807.0515l [hep-ph]]. 

P. Zemp, Pionic Hydrogen in QCD+QED: Decay width at NNLO, 
PhD thesis, University of Bern, 2004. 

V. E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 494 (2000) 9 
[arXiv:hep-ph/0009206| . 

G. Hdhler, in Landolt-Bdrnstein, 9b2, ed. H. Schopper (Springer, Berlin, 1983). 

J. Gasser, M. A. Ivanov, E. Lipartia, M. Mojzis and A. Rusetsky, Eur. Phys. J. 
C 26 (2002) 13 [arXiv:hep-ph/0206068] . 

U. G. Meissner, U. Ra ha and A. Rusetsky, Phys. Lett. B 639 (2006) 478 
[arXiv:nucl-th/0512035] . 

M. Hoferichter, B. Kubis and U. G. Meissner, Phys. Lett. B 678 (2009) 65 
[arXiv:0903.3890l [hep-ph]]. 

J. Gasser, M. E. Sainio and A. Svarc, Nucl. Phys. B 307 (1988) 779. 

A. Schmidt et at, Phys. Rev. Lett. 87 (2001) 232501 [arXiv:nucl-ex/0105010| . 

E. Korkmaz et ai, Phys. Rev. Lett. 83 (1999) 3609. 

C. Fernandez- Ramirez, A. M. Bernstein and T. W. Donnelly, arXiv:0907.3463 
[nucl-th] . 



[24] E. Fermi, Suppl. Nuovo Cimento 2 (1955) 58, K. M. Watson, Phys. Rev. 95 
(1954) 228. 



13 



